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During the period covered by this report. Professor Matzner did a detailed study of the 

/ 

effects of gravitational radiation on the relative positions of objects, 'with the hope that 
astrometric detection of gravitational radiation might be possible. His report is attached. 

The results are discouraging. It would appear that narrow-field instruments in jhe ATF 
class are still several orders of magnitude less accurate than would be required for this very 
delicate kind of measurement _ - r J — \ 

However, the situation changes considerably when wide-field instruments are considered. 
An instrument such as POINTS ought to be able to detect gravitational radiation at this level 
with ease. 

The list of bright quasars has been augmented. Of particular interest is the finding of 
quasars near two open clusters and one planetary nebula. These would be useful in 
determining absolute parallaxes of these objects. The report of the graduate stu dent who did 
this work is also attached. 
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GRAVITATIONAL WAVE DETECTION 
USING SPACE-BORNE OPTICAL INTERFEROMETERS 

The possibility of detection of long-period gravitational radiation using a high- 
accuracy (10 -5 arcsec), small-field (20 arcmin square) interferometer. 


The propagation of photons from astronomical and cosmological sources fol- 
lows null geodesics in the spacetime. Here we want to consider a flat space on 
which a propagating gravitational wave is superposed. The motion of a photon in 
the absence of the wave is described by a number of constants of its motion (es- 
sentially, the conserved components of its physical momentum). In the presence of 
the gravitational wave, there axe fewer, but still enough of, conserved momenta to 
solve for the photon orbit. However, the presence of the gravitational wave modifies 
the relationship of the conserved momenta to the physical direction of propagation. 
This can lead to deviation of pointing or to a displacement of the image in the 
image plane. 

Figure 1 shows a situation in which a net shift across the focal plane of the 
instrument will occur. The deflection of the target object is a function of the angle 
to the object referenced to the gravitational wave propagation direction, so objects 
separated by 66 in the image field will experience different deflections. However, 
the relative motion between two images is of order ah+(60) 2 , where h+ is the 
amplitude of the gravitational wave, a depends on the geometry and is typically 
small (see below), and 66 ^20 arcmin ^5 x 10 -3 (the field of view), so there is a 
strong suppression of the observability of the wave-induced deflection. 
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Figure 2 shows a possibility more appropriate to observing a real effect. 





Here, only the direction to the more distant source is deflected. The amount of 
deflection is then of order 0 / 1 +, where again o depends on the geometry and h+ is 
the amplitude of the wave. We now carry out the analysis and estimate angular 
deflection for various geometrical situations. 

For a plane wave travelling at the ^-direction, let 


u = z — t, 
v = z + t. 


Then, 


du dv = — dt 2 + dz 2 . 


The metric for weak gravitational waves is 

ds 2 = -du dv + [1 + h xx (u)] dx 2 + [1 - /in(tz)] dy 2 . (1) 

This is a typical h+ polarization pattern; the other (/ix^polarizafi 011 pattern in- 
volves a term h xy (u) dx dy that can be removed by a redefinition [a rotation in the 
(x-y)-plane]. The amplitude of the gravitational wave is h+ = h xx . 

In such a situation, there are still enough conserved momenta for photons to 
completely solve for the photon motion. (These are not the physical momenta, see 
below.) 

p x = const., p y = const., (2) 

Pv = +P<) = const - 
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If the gravitational wave were not present, p x and pt would be separately constant. 

The other component (not constant) of the momentum can be solved for be- 
cause the photon path is null 

9°^ Pa P/3 = 0 - 

To the accuracy required, using the inverse of the metric (1) expanded to first order 
in h, this is 

( Px + Pt)(Pz - Pt) + (1 “ h+)pl + ( 1 + h+)p 2 y = 0, (4) 


or 


(Pz-Pt) = 


(1 - + (1 + ft+)Py 


( 5 ) 


Pz+Pt 

The right side of this expression is constant except for the appearance of h. Using 
(3), we have, to lowest order in h+, 


b 1 
Pz = Pz + 2 


P\ -Pi 

b b 

Pz+Pt 


( 6 ) 


with the superscript “6” meaning “background” values (i.e., the values when /i+ = 0). 
Because of the form of the metric, we have p x = \p | cos 9 where 8 is the physical 
direction of propagation. (Here 9 is measured from the z-axis, the direction of 
propagation of the gravitaional wave.) Because of the wave, the physical direc- 
tion of propagation changes, since p x changes. From (6), we obtain the change in 
propagation angle 8\ 

A ^ h-|- sin 8 ^ cos 2^5 i 

&9 — ~r~ , /j i \') 

2 1 4- cos 9b 

b b 

where tan <f>b = Py/Px • 

There are also changes in the transverse ( <f > ) direction. Although p z and p y 
are separately constant even in the presence of the wave, they are not the phys- 
ical components. Instead, p x (physical) = 1/(1 + h + ) x ^ 2 p r and ^(physical) = 
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1/(1 - h+y/tpy, so that 


_ Py(physical) _ / 1 + h+ ( Py\ 

* " ^(physical) ]j 1 - h + y ? J 

« (1 + /i+) — = tan <f>b + /i+ tan <f>b, 


so that 


A^ = ft + (y 


tan <f>b 
+ tan 2 


As anticipated, both expressions for the deflection angle are proportional to h. 
In both cases, however, note that this deflection angle is a deflection in the photon 
path as it travels. Only if the wave is at the detector on Earth can this directly 
amount to a change in viewing direction. [This would be a subset of the case (a) 
considered above, and it yields diffferential deflection proportional to h+(A9) 2 , 
where A 6 is the angular offset between the two sources.] 


To analyze the apparent deflection in case (b) above, first note that the photon 
direction is the same after as before the wave passage. It is only during wave passage 
that the physical propagation changes. Hence, the situation can be idealized as a 
refraction in a plane sheet of glass, as follows. (jO/V't'l 

r w _ • 



Gravitational deflection angle 9^ 


Hence, as seen by the observer, the angular displacement of the source (compared 
to the no-wave case) is 

offset £9<^ 

A jT ___ Q — 

di + i cos 9 + d .2 d\ + c ?2 ’ 

Q 

assuming that the total distance to the source is large compared to length of the 
region occupied by the wave. 


4 


Since the angle ^is of order h, the resulting offset in the viewing angle 8 is of 
order 

(Atf)ob. 

where d is the distance to the optical source and £ is the duration of the gravitational 
wave pulse. In this expression we have dropped geometrical factors that axe typically 
of order unity. We have also inserted the average value of h+ throughout the wave; 
an oscillatory h+ will have a much smaller deflection. 

Because of the design of the space-borne optical interferometer, most sensitivity 
will be available for observed direction variations with period of order one year. A 
wave packet of ~ 1 light year scale in all directions (a single “positive-going” pulse 
with overall timescale ~ 1 year) is ideal. 

Ideal detection geometry would thus comprise two sources — one fiducial, one 
for detection — as close to the Earth as possible to minimize d. Candidates might 
be a pair of white dwarf stairs about ten light-years away and separated in distance 
by about one light-year. The suppression factor is then of order 1/10. 

Estimates of the amplitude h+ of the gravitational radiation can be made in 
several ways. A wave of period one year due to a binary star system (two 1/2 
solar mass stars, separated by 1 AU) would have amplitude h ~ 10 -8 /r(km) and 
so would be undetectable at any reasonable distance. A substantially relativistic 
collapse with a period of order one year would constitute a galaxy-mass collapse and 
is very improbable, but there is no clear single source candidate. Such a collapse 
would be detectable at 10 -n level from sources out to 300 Mpc. However, we can 
also imagine there is a stochastic background and estimate the amplitude of such a 
stochastic radiation field. A gravitational wave field has energy density proportional 
to (w/i+) 2 : 


p 

( w \ 2 

( h+ V 

Pclosure 

to 

vT 

ri 

V 10 -11 / 


Thus, if h+ ~ 10 -11 , then a uniform bath of random waves of period ~ 1 year would 
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be sufficient to “close” the universe. With the estimate of a ~ 10 -1 as the most 
favorable geometry, we would anticipate a maximum angular fluctuation ~ 10 -12 
radians, ~ 10 -6 arcseconds. This would seem to be about one order of magnitude 
below the design sensitivity of the instrument. 

The current observations of the millisecond pulsar give limits on the gravita- 
tional wave background in the ~ 1 year range. If we consider the recent report by 
J. Taylor as definitive, we find p/pdosure ~ 10 -4 (i.e., h+ ~ 10 -13 ), which gives 
at least another 10 2 suppression on the detectability below the capability of the 
instrument. We should note, however, that any single source (e.g., the millisecond 
pulsar) sample waves in only one particular direction. The interferometer proposed 
will have the ability to observe in essentially any direction. Furthermore, the grav- 
itational wave background can act as a noise source (diminishing with observations 
of more distant optical targets) most important for nearby stellar observations and 
should be considered in a signal- to-noise analysis. 
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More, dota search uoas eoodu cted during this semester . Quasars V-nousn 
as variables uoere included , and u)C took a closer look especially at: 2>c 

lOhieh is very tricot ard has a bright star nearby . 

W-t also studied -the list a| Hippare**’ Quasars and stars' . VTe picked the Wight 
orves and added them to our original list . V\)e limit the list to A5.*5 magnitude 
^j^sars . "they are all confined m "Table T . 1 , 


■from an object u>ith magnitude m ? photon ^lux can be calcoi «ted 
usirQ the -^ollocding refaTic?r\ 


S = N t X> z {a T") \o 


-0 4 m 


cohere : 


l\l*1fi 4 photons / (sec cm 1 nny) J^or a 0 magnitude AO star 
at a coav/e length ^ 5£o nm . 

t = transmittancy the top atmosphere totfhe detector 

C here ujt take to be unity ") 

Ik » diameter the telescope C t?5 x to 2- cm') 

Ad = bandpass the instrument used C dm on) . 

? 1 

Knowing thee repaired instrumental minimum e£ \o photons r toe c.an calculi st^ the 
ideal integirsttfv) time. The average total time coill he ©btamed by miitti plying the. 
ided integration time by Vito. Table "1 . a. e< plains. total observation t\me calcuA at" <fn . 


"the rej|t*rence stars are chosen coithin Ao arc minutes ^rjm eaelk quasar 
lOftk -the requirement that those objects: should bx, at Vast 5, arc sec . 

These re,|ertnc€ sWr are obtained ^rooi the Guide Star Catalog ^>ipace Spence 
XmtituV . K special computer program galled -pi elides d*ta then read ard VocaW them 
ground its quasar . 

Too o^en clusters and orvt planetary nebulae ushidh Arc uWich are ©lo^e. 
in pjst(<m to «ur quasars- or tlippaccrs- %W are Vsted in ""Talole It , 

M S throe der } ID-TJ., Astronomical Optics- / (Academic Thess Inc Sanlkfego , 10S tQ 


TABLE I Z 


Random Errors 

• Integration time needed to reach a given accuracy is increased by random errors. 



Error Sotirce 

El 

1. 

Photon Statistices (cro) 

1.000 

2. 

Background Light 

.002 

3. 

Image Shape/Size 

.30 

4. 

Image Motion (Jitter) 

(a) 

5. 

Grating Imperfections 

(a) 

6. 

Grating Motions 

(b) 

7. 

Grating Alignment 

TBD 

8. 

Field Modeling 

0.060 

9. 

Reduction Algorithm 

< 1 

10. 

Postfocal Response Variation 

TBD 

11. 

Reference Star Errors 

(c) 

12. 

Contamination 

TBD 


TOTAL TIME FACTOR: 

1.362 + TBD 


*F = Contribution to the integration time enhancement factor; 
total time enhancement is the sum of the individual F's. 

NOTES: (a) The design requirement corresponding to F 1 is feasible. 

(b) Included in jitter. 

(c) With proper selection of fields and reference stars, this error will be 
negligible ( F >C 1). 
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TABLE l 2(Cont.) 
Light Loss Effects 


• Integration time needed to reach a given accuracy is increased by light and 
other information losses. 


Source of Information Loss 

Type of Loss 

Throughput 

A. Grating Rejection 

Light 

0.25 

B. Mask for Grating Shadow 

Light 

0.75 

C. Grating Intrinsic 

Information 

0.50 

D. Loss in Optics 

Light 

0.50 

E. Detector Quantum Inefficiency 

Light 

0.10 

F. One- dimensioned Engine 

Information 

0.50 

G. Operational Interruptions 

Information 

0.50 


TOTAL THROUGHPUT: 

0.00117 


Integration time increased by 1 /throughput = 853. 

Observation Time Calculations 

• Total Observation = Integration Time for Ideal System 

X 

Time Factor for Light and Information Losses (853) 

x 

Time Factor for Random Errors (1.36) 

• For ATF: Overall Time Factor = 853 x 1.36 = 1160 

Observation Time = 1160x Ideal Integration Time 
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